VARIATIONAL EQUATIONS ON MIXED 
RIEMANNIAN-LORENTZIAN METRICS 



THOMAS H. OTWAY 



Abstract. A class of elliptic-hyperbolic equations is placed in the context 
of a geometric variational theory, in which the change of type is viewed as a 
change in the character of an underlying metric. A fundamental example of 
a metric which changes in this way is the extended projective disc, which is 
Riemannian at ordinary points, Lorentzian at ideal points, and singular on the 
absolute. Harmonic fields on such a metric can be interpreted as the hodograph 
image of extremal surfaces in Minkowski 3-space. This suggests an approach 
to generalized Plateau problems in 3-dimensional space-time via Hodge the- 
ory on the extended projective disc. Analogous variational problems arise on 
Riemannian- Lorentzian flow metrics in fiber bundles (twisted nonlinear Hodge 
equations), and on certain Riemannian- Lorentzian manifolds which occur in 
relativity and quantum cosmology. The examples surveyed come with natural 
gauge theories and Hodge dualities. This paper is mainly a review, but some 
technical extensions are proven. MSC2000: 35M10, 53A10, 83C80 
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1. Introduction: the projective disc 
In a small circle of paper, you shall see as it were an epitome of the whole world. 
Giambattista Delia Porta, 1589, on the camera obscura 

Analysis on Beltrami's projective disc model for hyperbolic space is in one sense 
very old mathematics. Beltrami introduced the projective disc in 1868 as one of the 
earliest Euclidean models for non-Euclidean space H] ; see also [TU] , |113| . But it also 
arises in the context of some new mathematics related to variational problems in 
Minkowski space and Hodge theory on pseudo- Riemannian manifolds. In practice 
Beltrami's construction amounts to equipping the unit disc centered at the origin 
of coordinates in with the distance function 

^ (l - y^) dx^ + 2xydxdy + (l - x^) dy^ 
(1 — a;^ — y^)^ 

Integrating ds along geodesic lines in polar coordinates, we find that the distance 
from any point in the interior of the unit disc to the boundary of the disc is infinite, 
so the unit circle becomes the absolute: the curve at projective infinity; c.f. |49j . 
Sec. 9.1. 

It is natural to ask how to interpret the so-called ideal points in the complement 
of the unit disc in K.^. It has been known for a long time that such points are not 
merely allowed by the projective disc model in a formal sense, but are actually 
useful in classical geometric constructions. For example, tangent lines to the unit 
disc can be used to characterize orthogonal lines within the disc, and certain families 
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of translated lines inside the disc attain their simplest representation as a rotation 
about an ideal point. Although these classical geometric operations on ideal points 
are well known [112j , [64] , geometric analysis on domains which include ideal points 
is still very incompletely understood. 

In this review we use Beltrami's model as a point of reference from which to 
survey aspects of geometric variational theory on mixed Riemannian-Lorentzian 
domains, on which the signature of the metric changes sign along a smooth hyper- 
surface. (The metric underlying ^ changes from Riemannian to Lorentzian on the 
Euclidean unit circle.) In Sec. 2 we consider variational problems which reduce to 
the existence of harmonic fields on the extended projected disc — that is, solutions 
to the Hodge equations 



on a domain of which includes the closed unit disc as a proper subset and is 
equipped with the Beltrami metric ([1]). Here d is the exterior derivative, with formal 
adjoint S, acting on a 1-form a. The crucial technical problem for the equations of 
Sec. 2 is the existence of solutions to boundary-value problems; this question is 
addressed in Sec. 3. Sections 4 and 5 are concerned with the related topics of 
duality and gauge invariance. These topics motivate the study of gauge-invariant, 
potentially elliptic-hyperbolic systems in fiber bundles. In Sec. 6 we briefiy review 
similar systems that arise in relativity and cosmology, and the relation of those 
systems to other elliptic-hyperbolic equations on manifolds which have been studied 
by mathematicians. 

On the extended projective disc the Hodge equations are no longer uniformly 
elliptic, but are elliptic on ordinary points inside the unit disc, hyperbolic on ideal 
points in the R^-complement of the disc, and parabolic on the unit circle, which is 
a singularity of the manifold. One might wonder why anyone would want to study 
such a peculiar system. There are at least two motivations for doing so: 

i) to learn what the geometry of the Beltrami disc reveals about elliptic- hyperbolic 
partial differential equations, and 

ii) to learn what elliptic- hyperbolic partial differential equations on the Beltrami 
disc reveal about the geometry of space-time. 

Regarding the first motivation, there is no canonical way to decide what consti- 
tutes a natural boundary-value problem for an equation that changes from elliptic 
to hyperbolic type on a smooth curve. Historically, physical analogies have been 
the main tool, chiefly analogies to the physics of compressible flow [T2]. However, 
it is also possible to approach the problem using a geometric analogy, which we will 
briefly describe. 

The highest-order terms of any linear second-order partial differential equation 
on a domain C can be written in the form 

(3) Lu = a {x, y) u^^ + 2/3(a;, y)u.j,y + "f{x, y)uyy, 

where (x, y) are coordinates on SI; a, (3, and 7 are given functions; u ~ u {x,y) . Tra- 
ditionally, the question of whether the equation is elliptic, hyperbolic, or parabolic 
has been identifled with the question of whether the discriminant 



is respectively positive, negative, or zero. If the discriminant is positive on part 
of fl and negative elsewhere on SI, then the equation associated with the operator 



(2) 



da = 5a = 



(4) 



A {x, y) = aj-' (3'^ 
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L is said to be of mixed elliptic-hyperbolic type. The curve on which the equation 
changes type is called the parabolic line or, borrowing the terminology of fluid 
dynamics, the sonic curve. A simple example of an elliptic-hyperbolic equation is 
the Lavrent'ev-Bitsadze equation 



for which the parabolic line is the x-axis. 

An alternative approach is to replace the class of differential operators L on a 
single domain of with the single Laplace- Beltrami operator, 



on a class of domains for which gij represents the local metric tensor. 

For example, the Lavrent'ev-Bitsadze equation can be associated to the Laplace- 
Beltrami operator on a metric which is Euclidean above the x-axis and Minkowskian 
below the x-axis. Elliptic operators will be associated to Riemannian metrics; wave 
operators will be associated to Lorentzian metrics. In this classification, the type 
of a linear second-order equation is not a function of the associated linear operator 
at all; that operator is always the Laplace-Beltrami operator. Rather, the type 
of the equation is a feature of the metric tensor on the underlying surface. Any 
change in the signature which results in a change in sign of the determinant g 
will change the Laplace-Beltrami operator on the metric from elliptic to hyperbolic 
type. The Laplace-Beltrami operator on surface metrics for which such a change 
occurs along a smooth curve will correspond to planar elliptic-hyperbolic operators 
in local coordinates. However, any curve on which the change of type occurs will 
necessarily represent a singularity of the metric tensor, as the determinant g will 
vanish along that curve. 

According to this point of view, in order to decide which boundary-value prob- 
lems are natural for a second-order linear elliptic- hyperbolic equation on R^, one 
should study the geometry of the underlying pseudo-Riemannian metric. For exam- 
ple, characteristic lines for the Hodge equations on extended can be interpreted 
as polar lines for a chord of the projective disc, and this suggests a natural class of 
boundary- value problems [95| . 

However, in this review we focus on the second motivation for studying harmonic 
fields on the extended projective disc. This is the opportunity that they afford to 
revisit, from a different point of view, work done by Chao-Hao Gu in the 1980's on 
the existence of extremal surfaces in Minkowski 3-space M^'-'^ [ll]-[l5], and similar 
variational problems which arise in the context of fiber bundles. 



The area functional for a smooth surface S in M^'^ having graph z = f {x, y) is 
given by 



The surface E is time- like when /J -f fy exceeds unity and space-like when /J -I- fy is 
exceeded by unity. Introducing Lagrange's notation p — fx, 1 = fy, the boundary 
between the space-like and time-like surfaces is the unit circle centered at the origin 
of coordinates in the pg-plane. 



sgn{y)uxx + Uyy = 0, 




2. The geometric variational problem 
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A necessary condition for S to be extremal on M^'-'^ is that its graph f {x,y) 
satisfy the minimal surface equation in the form [42j 

(5) (1 - qy + 2pqpy + (1 - g2) = Q. 

This is a quasilinear partial differential equation which is elliptic for space-like 
surfaces and hyperbolic for time-like surfaces. We can linearize this equation by 
the method of Legendre ([13], Sec. 1.6.1), applying the transformation z — px + 
qy — If {p, q) , X = ipp, y — ipq. We obtain the linear equation 

(1 - /) (fpp - 2pqifipq + (l - q^) fqq = 0. 

Adopting homogeneous coordinates {u,v,w) for w 7^ 0, we eventually obtain the 
equation [43j 

(6) [(1 - p') i'plp - 2p<zVp« + [(1 - q') y^,] q = 0, 

where p = —u/w and q — —v/w. Equation ([6]) can be interpreted as the Laplace- 
Beltrami equation on the extended projective disc P^; c.f. [!95^ and references 
therein. 

In models for more general classes of surfaces {e.g., Examples 2-4 from Sec. 2 of 
[109j ). the gradient of / may be replaced by the d-closed 1-form lu = ojidx + uj2dy. 
In particular, consider a space-like surface S embedded in M'^'^. The variational 
equations of the area functional can be written in the form 

(7) [(1 - w|) Wl]^ + {uJly + L02x) + [(l - t^l) ^^2] ^ = 0, 

(8) LU2x - ^ly = 0. 

These are the equations for an extremal surface which is space-like in some regions 
of M^^^ and time-like in others. However, the surface that it describes is singular 
on the circle + cj| = 1. 

If the Gaussian curvature of the surface is nonvanishing, then 

ljJlxl^2y — i^ly)'^ — ^lx^2y — ^ly^2x 7^ 

and the nonlinear system can be linearized by a Legendre transformation to obtain 

(9) [(1 - x^) u,]^ ~ {xyu,)y ~ {xyu2), -f [(l - y^) ^2]^ = 0, 

(10) Uiy - U2x = 0. 

This system has a geometric interpretation as the Hodge equations on the extended 
projective disc. 

From this point of view, solving the Hodge equations locally on extended is a 
way to approach generalized Plateau problems in Minkowski 3-space. For example, 
suppose one were to study the existence question for a class of area-extremizing 
surfaces in M^'-'^. The first step in such a program might be an existence theorem 
for a class of Dirichlet problems for eqs. (O, (fTU)) : The hypothesis of a sufficiently 
smooth vector-valued function h prescribed on all or part of a boundary curve in 
extended P^ would be shown to imply the existence of a 1-form satisfying the Hodge 
equations in the interior and equalling h on the boundary. 

Unfortunately, a good Hodge theory on extended P^ would not in itself provide 
a good variational theory even for stationary hypersurfaces in M^'^. One reason is 
that the Legendre transformation may itself introduce singularities. It is known 
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that, in the eUiptic region, such singularities can only occur on at most a point 
set if the image is the system dH), ([TOll [93] , However, higher-order singularities in 
the parabolic and hyperbolic regions of the equations are possible and, based on 
analogies to the equations of gas dynamics [22], expected. Moreover, although the 
Legendre transformation makes the equation simpler, it makes boundary conditions 
more complicated, so the interpretation of boundary-value problems for harmonic 
1-forms on extended P'^ in terms of extremal hypersurfaces in M^'^ is not generally 
straightforward. 

Nevertheless, the Dirichlet problem for a quasilinear elliptic-hyperbolic system, 
having a line singularity on the parabolic curve, is sufficiently formidable that the 
approach via linearization remains the one with the most apparent promise. (For 
a possible alternative, see |118j and references therein; for an alternative in the 
particular case studied by Gu, see [65].) 

3. Boundary- VALUE problems 

Both of the two motivations for studying the Hodge equations on extended 
begin with the existence question for boundary-value problems: under what condi- 
tions on the boundary can we expect solutions in the interior? The obstruction to 
answering this question is the problematic nature of boundary-value problems for 
elliptic-hyperbolic equations. 

The systematic study of partial differential equations which change from elliptic 
to hyperbolic type on a smooth hypersurface began in 1923 with the famous memoir 
of Tricomi [117] concerning the equation 

(11) yUa;x+Uyy=0. 

Tricomi's work was extended in the early 1930s by Cinquini-Cibrario 21j to include 
the equation 

(12) XUa:x + Uyy = 0. 

The study of these equations led to the establishment of two classes of normal forms, 
to which any linear elliptic-hyperbolic equation of second order can be reduced near 
a point (c./. [13j). and some appropriate boundary- value problems for such forms. 
In particular, the Tricomi problem for (|lip or (|12p prescribes the solution on the 
part of the boundary consisting of characteristic curves; see, e.g., Sec. 12.1 of [34] 
for a discussion. Tricomi problems for the scalar case of the system (jH]), (ITUl) were 
solved by Hua and his students in the 1980s [53], [55] , 

The boundary- value problems formulated by Tricomi and Cinquini-Cibrario were 
open in the sense that the boundary conditions are prescribed on a proper subset 
of the boundary. Closed problems, in which data are prescribed on the entire 
boundary, have only recently been studied in a systematic way [70] see also [103] . 
Closed problems were formulated as early as 1929 by Bateman [5], but in the 
informal style typical of British applied mathematics of that period. An isolated 
result, the existence of weak solutions to the closed Dirichlet problem for the scalar 
Tricomi equation, was proven by Morawetz in 1967 [81]; see also [104] . Some special 
cases are treated in [41], [71], [SH], and [115] . 
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The closed boundary- value problems in [70] are formulated for a class of equa- 
tions having the form 



where the type- change function K{y) is a continuously differentiable function sat- 
isfying certain technical properties, the most important of which are K{0) — and 
yK{y) > for y 7^ 0. In the simplest special case, K{y) = y, (|T3)) reduces to the 
Tricomi equation For this reason, equations of the form are said to be of 
Tricomi type. Similarly, eq. (|12p is the simplest example of an equation of Keldysh 
type — that is, an equation of the form [57j 

(14) K{X)UXX + Uyy = 0, 

where K{x) is generally taken to be a continuously differentiable function such that 
if (0) = and xK{x) > for x 0. Equations of this kind, with various lower- 
order terms, have arisen in the study of transonic gas dynamics and, recently, in 
optics; see, e.g., Sec. 3 of and Sec. 4 of [71], respectively. If the components of 
the vector (wi,M2) are sufficiently differentiable, then the system ((9]), (|10|) in polar 
coordinates can be represented as a single equation of Keldysh type. Sufficient 
conditions for the existence of unique solutions to closed boundary- value problems 
for equations having the general form (|14p are virtually unknown. 

The known results on the existence of solutions to boundary-value problems for 
dnj), pH]) can be summarized as follows: 

Let ri be the domain formed by the polar lines of a chord of the unit disc in 
extended P^. Then there exists a weak solution on fl with data prescribed on the 
non-characteristic part of the boundary. In fact, we can deform the chord in such a 
way that there is also a non-characteristic, explicitly hyperbolic boundary on which 
data are prescribed, providing a mild monotonicity condition is met. Solutions lie 
in a weighted function space. The weak solution is strong (in the sense of Sec. 
3.2, below) if we round off the sharp points on the boundary and perturb the 
lower-order terms. These assertions are stated precisely and proven in Sec. II of 
[95j . Although they are formulated in [95] in the context of projective geometry, a 
(different) boundary-value problem for an analogous domain in the scalar case has 
been formulated in the context of Minkowski geometry; compare Fig. 1 of |44| with 
Sec. 4, Fig. 2 of [Sg. 

However, we show in the following section that the closed Dirichlet problem for 
twice-continuously differentiable solutions of ([9]), pO|) is over-determined on the 
hyperbolic boundary. Note the odd terminology: Hyperbolic points are ordinary 
points lying inside the unit disc, and thus lie in the elliptic region of the eqs. ([9]). 
(|10|) . The hyperbolic boundary is the portion of the domain boundary on which the 
discriminant ^ of the equation is negative; so the hyperbolic boundary for ([5]), 
(|T0| lies outside the unit disc and consists of ideal points rather than hyperbolic 
points. 

3.1. The nonexistence of classical solutions. Expressed in polar coordinates, 
eq. ([6]) assumes the form 



(13) 



K{y)U^^ + Uyy = 0, 



(15) 
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for (j) = (j) {r,9) , provided r ^ 0. We define eq. (fT5|) over a sector 

= {(r, 0) |0 < ro < r < n, 01 < < ^2}, 

where 6*2 — 6*1 < tt and the interior of ilg includes a segment of the hne r — 1. For 
convenience we will take 9i to be negative and 02 to be positive. 
Writing eq. (llSp in the equivalent form 

(16) (1 - r^) 0„ + + r (1 - 2r2) 0, = 0, 

we show that the Dirichlet problem for this equation on a typical domain has 
a unique solution if data are given on only the non-characteristic portion of the 
boundary. This will imply that the Dirichlet problem is over-determined if data are 
prescribed on the entire boundary. 
Define the set £7+ C fig, where 

Q+ = {(r,6l) e R^lm < e < r < 1,-6*0 < 61 < 6I0} . 

We will choose the domain fl of eq. (fT6|) to be the region enclosed by the annular 
sector ^2+ and the intersecting hues tangent to the points (r, ±6) = (1, ±6*0) . 

Precisely, let {r,±9) = (l,±6'o) be the two points of intersection of the line 
X = xq, xq £ (e, 1) , with the unit circle centered at the origin of coordinates in R^. 
Let flo be the triangular region bounded by the vertical chord 70 given in Cartesian 
coordinates by 

70 = {ix,y) eR^\x = XQ,y^ < l-^;^} 
and the two polar lines of 70. (Recall that the polar lines of a chord are the tangent 
lines to the unit circle at its two points of intersection with the chord.) Then 

n = n+u no. 

In the following we denote by the subdomain of f2o consisting of ideal points, 
and by the arc of the unit circle lying between the points (r, ±6*) — (l,±6'o). 
Because e is a fixed number greater zero, the mapping from the region 17 in the 
polar r^-plane (the Cartesian xy-plane) to its image in the Cartesian r^-plane is 
well-defined, and we shall call the image H. as well. 

Theorem 1. Considering I116\) as an equation on the subdomain Q of the Carte- 
sian rO-plane, any solution cj) & (f2) of fl^i taking values f (r, 6) G (917) on 
the boundary segment 9f7\f7~ is unique. 

Proof. The method of proof has been applied to other elliptic-hyperbolic equa- 
tions [72], [I9|, [82]. Suppose that there are two solutions satisfying the boundary 
conditions. Subtraction yields a solution, which we also denote by (j), satisfying 
homogeneous boundary conditions. Define the functional 

rir.9) 

1=1 ipide + i)2dr, 

where 

(17) V'l = r2 (1 - r2) ct>l - (t>l 
and 



(18) i^2 = -I^Ae. 

on the Cartesian r0-plane. Because 
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(19) - 2(j)r [r^ (1 - r2) 0„ + (t>gg + r {l - 2r2) cj),] = 0, 

there is a function x (^i ^) such that xe — i^i ^^id Xr = fp2- 

The first step is to show that vanishes identicaUy on h'. Because (j) is the 
difference of two solutions having identical boundary values, (f){l,zL9o) = 0. Also, 
(pj. is zero along the lines 6 = ±9q. Then Xr = V'2 = ~'^'t'r4>e — on these lines. 
Integrating, x(''7~^o) = ci ^nd x(^i^o) — ^2, where Ci and C2 are constants. On 
the arc 

(20) xe - V^i = ~4>l < 

so ci > C2. On the intersection of the line r = e with dfl, 

Xe = ^1 = e (1 - e^) cj^j > 0, 

as ipe is zero there by the homogeneous boundary conditions, implying that C2 > ci. 
We conclude that ci = C2. Now (f20ll implies that = on t^. Integrating the 
differential equation xe = ~4'e ~ ^ using the boundary conditions (1, ±0o) = 
0, we conclude that vanishes identically on the parabolic line v. The homogeneous 
boundary conditions now imply that vanishes on all of dU,^ . The maximum 
principle for nonuniformly elliptic equations implies that <j) vanishes identically on 
ri+ (see the Remark following the proof of Theorem 3.1 in Sec. 3.1 of 37J). 

Denote by F the set of polar lines of 7, where 7 is the set of chords x ~ t for 
a:o < T < 1. Then F is a set of characteristic lines to eq. (fT6|) . which are intersecting 
pairs of tangent lines to the unit circle. Order F by decreasing radial distance, from 
the origin, of the point of intersection for the two polar lines of t. Then the indexed 
set Fa, beginning with the polar lines of 70 and having a limit point at the vertical 
line X = 1, foliates the subdomain in the sense that fi" can be expressed as the 
uncountable union of the elements of {Fa} . 

On elements of Fa we have the characteristic equation 

(21) dr^ - (r^ - l) dO^ = 0, 
so 

dx = XedO + Xrdr = (^xe ± Xrr\/r^ - l) dO 
and, using ^ and (fTSjl . 

= - (rv/r2 - IcPr ± < 0. 

Because Xe = Xr = on v, we know that x is constant on v. Integrating along 
elements of Fa and using the constancy of x on the arc v now implies that x is 
constant, and thus xe vanishes, on all of Q.^ . But xo = V'l ^nd 1 — < on 
so PT)) implies that 0^ = (/ig = on $7". Thus is constant on fi^. Because 
vanishes on we find by continuity that the constant is zero. 
This completes the proof of Theorem 1. 

Corollary 2. Considering \lb]) as an equation on a domain VL of the Cartesian 
rO-plane, the closed Dirichlet problem for classical solutions of eq. il6]) is over- 
determined on fl. 
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Remarks, i) A nonexistence result of an entirely different kind, for solutions of 
a semilinear elliptic-hyperbolic equation, is presented in |68j . That result is based 
on an elliptic-hyperbolic extension of the Pohozaev identities of elliptic theory. 

ii) Defining the first-order operator N(j) = —20^ and the 1-form 

eq. (|19p can be interpreted formally as associating a conservation law to the second- 
order operator L(f) defined by eq. (|f 6p via the method of multipliers: 



= J J N<j)- L(j) drde = J J dip. 

(Compare this equation with p. 262 of [69] and eqs. (5), (6) of [80]; see also the 
equation preceding eq. (7.1) of [29], taking dU/dt = in that equation.) However, 
this interpretation is only applicable to a sufficiently small coordinate patch about 
an arc of the sonic curve. In general, x is only incompletely specified by eq. (|19p : 
however, in the case considered in Theorem 1, regions on which (j) must vanish are 
deduced from boundary conditions and the particular geometry of the domain. 

in) Despite the close similarity of Theorem 1 to Sec. 3 of [82], eq. (fTH]) differs 
from the corresponding equation, (23), of [82j in its type-change function and lower- 
order terms. Moreover, because the hyperbolic boundary depends on the form of the 
characteristic equation, this difference in type-change function affects the geometry 
of the domain. 

iv) The characteristic equation (|2ip . in its Cartesian form 

(22) (a^ - x^) dy^ + 2xy dx dy + (a^ - y^) dx"^ = 

for a a constant, seems to have been introduced in 1854 by G. G. Stokes, in the same 
Cambridge Smith's Prize examination in which Stokes' Theorem first appeared 
(c./. [114] and eq. 1 of ^S]). Although no geometric context was suggested in the 
examination question. Stokes drew the attention of examinees to the possibility of 
singular solutions. 

3.2. A strongly well-posed boundary- value problem. Solutions which exist 
in the closure of the graph of the differential operator are said to be strong; c.f. [33j , 
p. 354. Precisely, a strong solution of a boundary-value problem for an operator 
equation Lu = /, with / G L^, is an element u G for which there exists a sequence 
u'^ of continuously differentiable functions, satisfying the boundary conditions, for 
which 

lim \\u^ — u\\t2 = lim \\Lu^ — /Ilr2 — 0. 

v — ^oo u — >oo 

A consequence of this definition is the uniqueness of strong solutions. 

A system of two differential equations on a domain of is symmetric positive 
[33j if it can be written in the form of the matrix equation 

(23) Lw = A'^Wx + A^Wy + Bw = Lp, 

where w — {wi,W2) is a vector, the matrices A^ and A^ are symmetric, and the 
symmetric part k* of the matrix k = B — (1/2) {A], + A^^ is positive definite. Here 
and below, an asterisk denotes symmetrization: n* = {k-\- k^)/2. 

Considering a solution to to (|23p as a mapping from the domain to a finite- 
dimensional vector space V, denote by dV the restriction of V to the boundary dU, 
of f2. We will use in a fundamental way a theorem by Friedrichs |33j : see also |63j 
and [1(17] . 
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Friedrichs' Theorem on Symmetric Positive Systems. Let fl be a bounded 
domain o/M^ having boundary dfl. Consider the matrix j3 — niA}^Q^ + n2A^Q^, 
where n — (ni,n2) is the outward-pointing normal vector to dfl. Suppose that the 
matrix {3 admits a decomposition (3 = (3++ j3-^ for which the sum of the null spaces 
for 13+ and (3^ spans dV\ [H+ n$H_ = 0, where are the ranges of I3±; the matrix 
l-i ~ f3-^ — [3- satisfies fi* > 0. Then the boundary-value problem given by 

(24) Lw ^ ip in n, 

j3-w = on (9il, 

where L is given by I123\} . has a strong solution whenever the system is symmetric 
positive and the components of if are square- integrable. 

A boundary- value problem which possesses a strong solution is said to be strongly 
well-posed. We show the existence of such a problem for a class of equations of 
Keldysh type. We initially consider equations having the form 

(25) [K{rj)u,l^ + u^^ + ku^ = f{i^,i), 

where K'{r]) > 0, A; is a nonzero constant, if (77) < for < 77 < rjcrit and if (77) > 
for rjcrit < 77 < i?- We do not expect classical solutions, so rather than study this 
equation directly, we consider the associated system 

(26) Lw = £-Wn + A^W(^ + Bw = F, 

where L is a first-order operator, w = (^1(77, ^), k;2(77, 0) i -f" — (/7O) i 

. ( » ) , ^ ( ; ; 

and 

(28) (*■'•'" I 

We interpret 77 as a radial coordinate and ^ as an angular coordinate, so that the 
system (P5| - (|28p is defined on a closed disc of radius R. The system is equivalent 
to (j25[) if the components of w are and wi — u^i, W2 — u^. 

Theorem 3. Suppose that there is a positive constant i^q such that K' [rj) > vq. 
Let there be continuous functions (t(^) and t(^) such that the boundary condition 

(29) (7(0^1+^(0^2=0, 

where the product cr(0''"(0 *s either strictly positive or strictly negative and has 
sign opposite to the sign of k. Then the boundary-value problem Il26\) - i29\) . with 
K and k as defined in eq. \25\) . possesses a strong solution on the closed disc 
{(^jO |0 !i ?7 ^ i?} provided \K{0)\ is sufficiently small. 



Proof. Multiply the terms of eq. (p6|) by the matrix 



c a 

where a and c are constants; the sign of c is chosen so that arc < (so cfc > 0); 
a > 0; and |c| is large. The matrix E is nonsingular provided 

det£: = a'^ + c^K(7j) ^ 0. 
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Because K{r]) is continuous, increasing, and K^rjcrit) — 0, the invertibility condition 
for E becomes 

min \K(rj)\ < —. 

This condition wiU be satisfied provided |i^(0)| is sufficiently small. 
The symmetric part k* of the matrix 

K^EB- {1/2) [{EA\ + {EA^)^] 

has determinant 



ak 
~2 



cK'{r]) 



ak 
Y 



> 



ak 



ak 
T 



so the resulting system is symmetric positive provided |c| is sufficiently large. The 
proof will be complete once we show that the boundary conditions are admissible. 
Proceeding as in [116j , choose the outward-pointing normal vector to have the form 
n — K^^(Ti)dr\. 

Then on the boundary rj = R, 



f3 = 



c 



(i?) 



Choose 



1 



Choose /3+ = (3 — (3-. Then (3-w = 0, as 
circle rj — R. Moreover, 

M = 

1 / (t^ - a^) a-2(7Tc 

2aT\aK-^{R)\ 



implying that 



(72) C- 



T c + ma 
~T^aK-^{R) +aTC 

implies that W2 — — (cr/r)wi on the 



((T^ — r^) c — 2ara 
- a^) aK-\R) - 2aTC 



M = 

(r2 - cr2) a - 2<7TC ara {K-^{R) - l) 

cr-+T- y aTa{K-^{R) -1) {t^ ^ a^) aK-\R) - 2aTc 
If (TT < 0, choose c > 0; if err > 0, choose c < 0. Then the matrix will be 
non-negative provided |c| is sufficiently large. 
Now 



aUli + TUJ2 



and 



TWi — aUJ2 
0-2 -I- t2 



TC -f- era 

-TaK-^{R) +ac 



Ta — ac 
aaK-^{R) + TC 



so dK- n d\+ = 0. Because conditions are given on the entire boundary of the disk, 
the null space of /3_ alone spans the range of dV. 

The invertibility of E completes the proof of Theorem 3. 



Remarks, i) Equation ([25)1 is a lower-order perturbation of an equation stud- 
ied by Magnanini and Talenti in the context of singular optics [7T]. Those authors 
considered data prescribed on the boundary of the unit disc. They were able to 
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construct an explicit, smooth solution in the interior of the disc having a point sin- 
gularity at the origin. Moreover, they showed the problem to be weakly well-posed 
on the entire domain. Theorem 3 of this section replaces Theorem 3 of j95j which 
claims, on the basis of an erroneous proof 96J, that a strong solution to a closed 
boundary-value problem exists for an arbitrarily small lower-order perturbation of 
the equation studied in '71] . In the present theorem the perturbation is lower-order 
but not arbitrarily small. 

ii) A similar result has been proven by Torre [116j for the equation 

(30) ^{vu,\+(^^-^'y^^^f{v,0, 

where / is a sufficiently smooth function and w is a constant. In [116] . condition 
((29| is imposed on functions a, r such that the product ctt does not vanish on 
the outer boundary of an annulus and cr, r satisfy the Dirichlet conditions ti = 1, 
T = on the inner boundary. Because eq. I|30p is a helically reduced wave equation 
in M^'^, it provides a toy model for the reduction of the Einstein equations by a 
helical Killing field (c./. Examples 1 and 4 of Sec. 6, below). In distinction to eq. 
(l25l) . which is of Keldysh type, eq. ((30)) is of Tricomi type. 



This section has been devoted to closed boundary-value problems for harmonic 
fields on extended and similar equations of Keldysh type. The vast literature 
on open boundary- value problems for equations of the form (jl3p and its nonlinear 
relatives has been ignored. The literature on such problems published during the 
first half of the twentieth century is reviewed in [T2J; for more recent literature, see 
[T7] . See also [72] . All three references concentrate on the literature in mathematical 
fluid dynamics, but results also exist on boundary-value problems for nonlinear 
equations related to (|13p that arise in connection with the isometric embedding of 
Riemannian manifolds into a higher-dimensional Euclidean space [67], [84], |127j . 

4. Duality 

In this and the following section we return to the geometric variational problems 
introduced in Section 2. But we will eventually replace the linear Hodge theory of 
eqs. (12) by a nonlinear Hodge theory based on an elliptic-hyperbolic extension of 
the Yang-Mills equations. 

4.1. Hodge duality of surfaces. It has been known for nearly a century that 
every solution of the Euclidean minimal surface equation 

(31) V-f , ^" Uo 

over all of is an afhne linear function fVP . Geometrically, this means that a non- 
parametric minimal surface of M.^ is a plane. This Bernstein property was extended 
to entire solutions of the Lorentzian maximal space-like hypersurface equation 

(32) V-f , ^" Uo 

by Calabi [H]: a maximal space-like hypersurface of M^'^ is a space-like plane. 
Calabi's result was later extended to space-like hypersurfaces of n-dimensional 
Minkowski space by Cheng and Yau [20] . 
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Yang has interpreted this kind of equivalence between equations ([3T|) and ([32]) in 
the language of differential forms ( [126j . eqs. (2.23)— (2.29); see also [1]). In Yang's 
interpretation, Calabi's equivalence becomes duality under the Hodge isomorphism. 

Yang's argument involves writing eq. ([5]) in the dual form [108) 

(33) d [p{Q)uj] = = 0, 

where, as in ([21), d : A^" ^ Ap+^ is the exterior derivative with formal adjoint 
(5 : AP -> AP-i; 

* : A^ ^ A"~P is the Hodge involution. If the surface is embedded in Euclidean 
3-space, then we choose 

(34) piQ) ^ 



VT+Q 

If the surface is embedded in Minkowski 3-space, then 
(35) p{Q) = ^ 



In the latter case it is necessary to distinguish space-like surfaces, for which Q < 1, 
from time- like surfaces, for which Q > 1. The surface Q = 1 is the light cone. Note 
that none of these surfaces can be found in M^'^ without solving the system 

If a; is a p-form satisfying (|33p over a region having trivial de Rham cohomology, 
then the n — p — 1-form * [p{Q)duj] is closed, and thus is the differential of an 
n — p — 2-form a : 

da = * [p{Q)diLj] . 
If p{Q) is given by p4)) . then this implies that 



= d^LV 



(-1) 



da 



* 



v/1 - \d<y\ 



±S [p{Q)da] , 



with p(Q) given by see eqs. (2.23-2.30) of [T26j for details. 

Yang obtains from this duality, not the Bernstein property, but rather a Liouville 
property ( |126] . eq. (2.31)): if either dui or da have finite energy on K" in the sense 
that 



(36) E = / p{s)ds * 1 < oo 

JR" Jo 

for p{Q) given by either pi]) or ([35| . then 

(37) duj = or da = 0. 

See Sees. 1-4 of [llOj for more details of this and similar Liouville theorems. We 
return to this topic in Sec. 5. 

4.2. Hodge duality of flow metrics. It is also possible to establish a kind of 
Hodge duality between extremal surfaces in M^'^ and a fully quasilinear model for 
compressible potential flow. If L is the operator of eq. ([3]), we define as in [T^] the 
flow metric for the operator L to be the metric tensor gij having distance element 

ds^ = a (x, y) dy^ — 2/3 (x, y) dxdy -t- 7 (x, y) dx^ . 
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The metric gij is Riemannian for regions on which L is an elhptic operator and 
Lorentzian for regions on which i is a hyperbohc operator. 

This idea does not use the hnearity of L in any essential way, and its original 
motivation seems to have been the quasilinear equation for the steady flow of an 
ideal gas. Define the quantity 

where u and v are the components of the velocity of an ideal compressible fluid in 
the X and y directions, respectively; 7 > 1 is the adiabatic constant of the medium. 
(The notation, which is historical, implies that the flow metric has traditionally not 
been applied beyond the cavitation velocity u'^+v^ =2/ (7 — 1) . The gas is assumed 
to be adiabatic and isentropic in order to avoid having to specify thermodynamic 
variables.) 

In the irrotational case we can deflne a potential function ip (x, y) for which 

dij} = udx + vdy 

and 

— udy — vdx^ 

where the centered asterisk again denotes Hodge involution. Then the flow metric 
for the continuity equation of compressible ideal flow in the plane, 

— U^) Ux — 2UVUy + (C^ — U^) Vy — 0, 

has distance element 

(38) ds^ = {dx^ + dy^) - {*dijf ; 
the flow metric of the classical minimal surface equation 

(39) (1 + w^) - 2uvuy + (1 + u"^) Vy ^0 

for the graph of a smooth function z = (x,y) has distance element 

(40) ds'^ ^ dx^ + dy^ + diP^ . 

Each of the two flow metrics (|38p and (1^0]) decomposes into a conformally Eu- 
clidean part and a non-Euclidean part. Moreover, if the classical minimal surface 
equation ((39)) for smooth surfaces embedded in is replaced by the equation ([5]) 
for surfaces embedded in M^'^, then the non-Euclidean part of that flow metric is 
the Hodge dual of the non-Euclidean part of the flow metric for ideal compressible 
fluids: the flow metric for (O has element of distance ds" given by 

ds"^ = dx^ +dy^ - dip^. 

For more details on duality of flow metrics, see Sec. 2.1 of [91] . 

The approach to duality initiated by eq. ([55]) seems to have been introduced in 
[108] . in the context of compressible flow on a compact Riemannian manifold; also 
compare [51] and [75] with [TUS]- 
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5. Gauge invariance 



5.1. Twisted nonlinear Hodge equations. It is natural to ask why the Bern- 
stein property for extremal surfaces in [TT] , [TS] , and [5D] becomes a Liouville prop- 
erty for differential forms in eq. (|37p . The Bernstein property for solutions of the 
two equations pip and (|32p has an intuitive interpretation because these partial dif- 
ferential equations can be embedded in a geometric theory, the theory of extremal 
hypersurfaces. In order to obtain an intuitive interpretation of Yang's Liouville 
property for solutions of the two systems , ([51)1 and , (|55|) , it would be use- 
ful to embed those partial differential equations as well in an appropriate geometric 
theory. Yang observes |126j that we can do so if the order of the differential form 

is 1 or n — 1. We can also do so in the case studied by Yang, in which the order 
of the differential form is 2. 

In that case eq. with p given by ([55]) . is the Born-Infeld equation for an 

electromagnetic field lu [T3], [126] . There is a generalization of the Born-Infeld 
equations which illuminates their geometric properties and, in particular, gives a 
geometric interpretation of the Liouville property for those equations. 

Proceeding as in [89 and 90J , we denote by X a vector bundle having compact 
structure group G C SO{m) and define X over a smooth, finite, oriented, n- 
dimensional Riemannian manifold M. Form an admissible class of connections by 
choosing a smooth base connection D in the space of connections compatible with G 
and considering the class of connections D + A, where A is a section of adX®T*M 
which lies in the largest Sobolev space for which the energy functional E given by 
(j36| is finite (with the domain R" replaced by Af); details are given in |119j for the 
case p = 1. Denote by Fa the curvature 2-form associated to a connection 1-form 
A. Then Q = \Fa\'^ — {Fa, Fa) is an inner product on the fibers of the bundle 
adX (g) (T *M) . The inner product on adX is induced by the normalized trace 
inner product on SO{m) and the inner product on (T * M) is induced by the 
exterior product * {Fa A *Fa) ■ Sections of the automorphism bundle AutX are 
gauge transformations. These maps act tensorially on Fa but affinely on A; see, 
e.g., Sec. 6.4 of [75]. Note that it is possible to describe the geometry of Fa either 
in terms of a principal bundle or in terms of the associated vector bundle. The 
choice is a matter of convenience; see, e.g., Sec. 2.3 of T3J for a discussion. 

In the smooth case we take variations by computing {d/ dt){F£)-^-tA) at the origin 
of t. Using the fact that for any smooth section a we have (^SSj, Sec. 1) 



FD+tA{(T) = {F + tDA + i 



'AAA) (a). 



we obtain 




(41) 




Letting t — 0, the right-hand side of (PT|) can be written 



(42) 
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We assume that either dM = or, if not, that Fa satisfies a "Neumann" boundary 
condition of the form 

(43) r(*F)=0 

on dM, where i* is the pull-back under inclusion of the boundary of M in M. 
This is equivalent in local coordinates to prescribing zero boundary data for F in 
a direction normal to dM; see [74] for details in the case p=l. 
Set SE equal to zero. Then (g!]) and (gl]) imply 

= / {DA, p{Q)Fa) dM = 

J M 

f A^ A{p{Q)Fa)j, + f {A, D* {p{Q)Fa)) dM, 

JdM JM 

where D* denotes the formal adjoint of the exterior covariant derivative D; -d de- 
notes tangential component on the boundary and N , the normal component there. 
Condition (|43|) implies Euler-Lagrange equations of the form [87] , [88] 

(44) D* {p{Q)F) = 0. 

Because F is a curvature 2-form, it satisfies an additional condition 

(45) DF = 0, 

the second Bianchi identity. 

li p = 1, then eqs. (|44p . (|l5)) degenerate to the Euclidean Yang-Mills equations 
of high-energy physics. Viewed in another way, one obtains eqs. (g^ . from 
eqs. ((33l) by twisting the cotangent bundle in which solutions of (|33)) live. This 
twist is represented analytically by a nonvanishing Lie bracket with A. Precisely, 
A £T (M, adX ® T*M) is a connection 1-form on X having curvature 2-form 

Fa ^ dA+]^[A, A] ^ dA + A h A, 

where [ , ] is the bracket of the Lie algebra 3, the fiber of the adjoint bundle 
adX. Thus eqs. (j44|) . (|45|) can be characterized algebraically as a non-commutative 
version of p3p , or geometrically as twisted nonlinear Hodge equations. 

Equations (|44|) . ([45]) change from elliptic to hyperbolic type along the curve 

So, in this case as well, the elliptic and hyperbolic regions of the equations cannot 
be found without solving the system. For p given by eq. ((34|) . eqs. ((44|) . ((45|) do not 
change type, but ellipticity degenerates as Q tends to infinity. 

Because eqs. are equivalent to eqs. (gl)) . (g5|) in the special case of an 

abelian gauge group G, the closed 2-form lo of the classical Born-Infeld theory is, 
geometrically, an abehan special case of the curvature 2-form Fa [c.f. [90], Sec. 
1). We can interpret Yang's Liouville theorem given by eq. (IST]) for solutions of 
the system ([33]) with (p4]) or ([SSjl . in the generalized perspective of [88]- [90], as 
the assertion that a finite-energy solution of ([55]) on R" with density or ([55]) is 
associated with a bundle having zero curvature, just as entire solutions in M?" or M^'^ 
are associated with surfaces having zero curvature. That geometric interpretation 
is reflected in the degeneration of the Bernstein property of [15j , for entire solutions 
of ((3T]) and ([32]) . to a Liouville property in [126j for entire solutions of ((33]) with 
dSl) or ([35]). 
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Liouville theorems have also been derived for solutions of (j44| and similar geo- 
metric objects; see, e.g., [86], [S^, Sec. 5 of [93], and Sees. 3 and 4 of |110j . An 
n-dimensional Born-Infeld system related to time-like extremal surfaces in higher- 
dimensional Minkowski space is studied in |59| . 

In [56| a physical model for compressible ideal flow with 50(3) gauge invariance 
is proposed. The mathematical structure of this model is similar in some ways to 
the variational theory presented here. 

Choosing the function p{Q) in eqs. (|^^ to be 1 -I- AjQ^-^, where 7 is a param- 
eter, the gauge group to be the abelian group C/(l), and the underlying metric to 
be Robertson- Walker, eqs. (|44|) become identical to eqs. (12) of a nonlinear 
electrodynamic model for the accelerated expansion of the universe. 

A version of eqs. (|44|) . (f45|) for mappings exists [9^, [94], [9^ which is related to 
harmonic maps in exactly the same way that eqs. (|44|) . (|45|) are related to Yang- 
Mills fields. This class of maps is a special case of the F -harmonic maps introduced 
by Ara [5]. 

5.2. Global triviality on R". In this section we extend the Liouville theorems 
cited in the preceding section to the case of an energy functional for which the 
variational equations may be elliptic-hyperbolic and the density may cavitate, in 
a pointwise sense, on a set of measure zero. The bundle connection need not be a 
pointwise solution of any equation at all in order to satisfy the hypotheses of this 
section. 

Consider a principal bundle 11 over a domain fl of M" and a 1-parameter family 
(/)* of compactly supported diffeomorphisms of H. for which o 0* = with (fP 
the identity transformation. This family of reparametrizations can be lifted to the 
principal bundle by parallel transport, with respect to an arbitrary smooth connec- 
tion, along the curve Xg = (f'^ix) from xq = x to Xt = 4>^{x). li A gT (fi, A^ (adll)) 
is a Lie-algebra- valued connection 1-form, then we define A* = (V'*)* A, where ip* is 
a lifting of (p* to F (fi, A^ (adll)) . (The superscripted asterisk denotes an induced 
mapping.) For details of this lifting, see |106j . We say that the connection A is 
r-stationary with respect to the energy functional E given by (l36|) if 



c.f. [2j. The following theorem provides sufficient conditions under which A must 
have zero curvature almost everywhere if $7 = E". A slightly more complicated 
argument will extend the result to manifolds of constant negative curvature {c.f. 
[106] ) and to ]R"\E, where E is a singular set of prescribed codimension {c.f. [5S]). 

Theorem 4. Consider the geometric construction of the preceding paragraph. 
Assume that p is continuously differentiable, nonnegative, and weightless under 
conformal transformations; that there is a positive number Qcrit such that 





(46) 
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with Q(F) ~ \Fa\'^ defined as in Sec. 5.1; that the restriction of E to a Euclidean 
n-disc of radius R is r-stationary Vi? > 0; that p'{s) < OVs > 0; and that 

(47) E^B^ < CR'^ 

as R tends to infinity for a positive constant C and a sufficiently small (positive) 
constant k. Then if Q is hounded above by Qcrit, Q iF{x)) must be zero for almost 
every x e M" for n > 4. 

Proof. Similar theorems under somewhat different hypotheses have been proven 
in |106j : see also Sees. 2 and 5 of [SS], and Sec. 5 of [93]. Assume that the lifting 
described earlier has been constructed. Denote by V'* the lift of a 1-parameter 
family of compactly supported diffeomorphisms of Q such that 

and 'ip'^ = identity. Define 

(48) f = i;\x)^x + tax)+0{t^), 
where 

at |t=o 

is the variation vector field, to be chosen. We have 



^ F^Af)d^df^ = ^ F,Af)^dx''^dx'" 
dt\t=o ' dt\t=o dxf" 



dp 



dx^' 



dt\t=o 



F^Af) [ siSL + sUg^ + SLt^ + o{t^) ) dx^dx 



dx'' 



Make the coordinate transformation a; — > y, where 

y={^'y'{x). 

Then 



A^^F.,(/)d/M/^=2F.,(.)|; 



dx^dx^ . 



If J is the Jacobian of the transformation a; — > y, then 



at |t=o 

By hypothesis, Vi? > 1, 



d 

dt \t=o 



dx 



div ^. 



dt \t=o 



e{Q)- J 
dt |t=o 



df 



* 1 



* 1+ 



(50) 



e'{Q)2 l^j^^^^FAf)df'df=,Firn{f)df'df'^^ * 1 
Substitution of dH]) into ^ yields 

e{Q)divi* 1 = 4^ e'{Q) (^F,,{x)^dx'dx^ , F,,{f)df 



'dp ) * 1. 
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_d_ 



Choose an orthonormal basis 

r _ \ d d 

and let |106j ^ = r](r)r ■ d/dr, where: r is the radial coordinate in a curvilinear 
system; ri{r) G [0, 1] ; r]'{r) < 0; r]{r) — v{r/T) = 1 for r < t, where r is 
a number in the interval (0, 1); there is a positive number S for which r]{r) — 
whenever r exceeds t + S. For this choice of ^, 

2 



(51) 



e{Q) (nrj + rrj) * 1 



Qp{Q)n * 1 



* 1. 



Our hypothesis on the sign of p' implies that 

^« d , 



Qp{Q) 



ds 



sp{s)) ds 



(52) 



[sp'{s) + p{s)] ds < I p{s)ds = e{Q). 



Substitution of ^ into ^ yields 

e{Q) {n-q ~ 4r] + rrj') * 1 < 4 

By construction 



p{Q)rij' 



* 1. 



rrj' [r) 



This yields 



< 4 



B, 



d_ 

dr 



d_ 



* 1 < 



< 0. 



e(Q) 



Br 



(4 — n)i] + 'T-Q^^v 



* 1. 



As 5 tends to zero we obtain 

< (4-71 



Multiply this last inequality by the integrating factor t'^ and integrate over 

T between ri and r2. We find that 

„4— n ; 



(53) 



We can write the growth condition (I47p in the form 



,4+/c-n 



(54) 

where 4+A; — n < for sufficiently small k. The right-hand side of (|54p tends to zero 
as r tends to infinity. The left-hand side is nonnegative by construction. Thus the 
conformal energy r'*~"i?|B^ tends to zero on R". Because by ([55)) the conformal 
energy is nondecreasing for increasing r, we conclude that E is identically zero on 
K". The vanishing of the energy on a ball of infinite radius implies the pointwise 
vanishing of Q almost everywhere by inequality (|46p . This completes the proof. 



An exactly analogous result holds in the case of a velocity field for a steady 
polytropic flow of an ideal compressible fluid; see Sec. 5 of [S3] for details. 
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5.3. Local existence and regularity. In this section we briefly review the fun- 
damental analysis of eqs. ((44|) . (f45|) . For details, see [89]; see also Sec. 2 of [93] . 

The local existence of solutions has been shown in the case p{Q) = Q*^^"^^/^ 
for p exceeding n/2, where n is the dimension of the domain. In this case the 
energy functional is Palais-Smale, and the existence of weak solutions follows by 
conventional variational theory. The equations are nonuniformly elliptic for this 
choice of p. It is possible to show that, in general, any weak solution is Holder 
continuous on compact subdomains whenever Fa G for p > n/2 and p is chosen 
in such a way that the system is uniformly elliptic. 

In the special case of an abelian gauge group, a weak solution of ()44[) . (|45|) is any 
curvature 2-form Fa for which p{Q)Fa is orthogonal in to the space of d-closed 
2-forms dC, € L^{B) such that C € has vanishing tangential data on dB. For a 
nonabelian gauge group, an obvious extension to inhomogeneous equations allows 
a weak solution to be defined by the equation 



Here i? is a Euclidean n-disc, so the natural hypothesis on the domain is that it 
can be covered by n-discs; for example, assume that the domain boundary has no 
cusps. 

We briefly outline the proof of regularity: The first step is to derive a weak 
subelliptic estimate for using difference quotients to establish the existence 

of an _ff ^'^-derivative. The subelliptic estimate allows one to conclude from a lim- 
iting argument that |F| is locally bounded by Morrey's Theorem. Now we apply 
a mean-value inequality for Lie-algebra-valued sections, which is valid in the uni- 
formly elliptic case of the equations (|44p. (|45|) . The mean- value inequality is applied 
to points in the solution space, so that the L^-diflference of weak solutions can be 
compared despite the nonlinearity of p. Using this inequality, it is possible to esti- 
mate the L^-difFerence between a 2-form weakly satisfying eqs. ((33)) and a bounded, 
weak solution of eqs. (04]), (|45]) . The latter is considered in an exponential gauge, 
fixed in a Euclidean n-disc B centered at the origin of coordinates in R". Solutions 
to eqs. p3p have known regularity, and in an exponential gauge, A[Q) = and 



Thus elliptic estimates, and the minimizing property of the variance by the mean 
with respect to location parameters, eventually show that the i^-difference of F 
and its local mean value decays sufficiently rapidly as the radius of the underlying 
domain shrinks to zero that Holder continuity can be obtained from the Campanato 
Theorem. The final step is to show that the Campanato estimate is preserved under 
continuous gauge transformations in a small n-disc centered at a point close to the 
origin. This allows one to apply a covering argument which will extend to the entire 
domain the estimate that was obtained in an exponential gauge at the origin. 

5.4. Other symmetry groups. Although the system (I44p . (05]) is invariant under 
the action of the gauge group G and has the potential to change from elliptic to 
hyperbolic type, the system is so poorly understood that it is hard to know how 
these two properties affect solutions. In principle, the system cannot even be said 
to be elliptic- hyperbolic, as the type of the equation can only be fixed modulo 



(55) 




\A{x)\<-\x\- sup |F(y)|. 
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gauge transformations. However, the gauge invariance can be broken and ()44|) . 
()45p represented as an elliptic-hyperbolic system whenever the curvature Fa lies 
in a sufficiently high Sobolev space. It is this property, which is shared with the 
Yang-Mills equations [119j . that allows one to formulate sufficient conditions for 
the existence of solutions, and their regularity in the uniformly "subsonic" case, as 
described in the preceding section. 

The distance element on the Beltrami disc is also invariant under the action of 
a gauge group, the group of projective transformations. This invariance has the 
consequence that the parabolic line for the Hodge equations ([5]), pU|) on extended 
is gauge-invariant under perspcction to any conic section. So ([9]), (jlO[) is an example 
of an elliptic-hyperbolic system in which gauge invariance has a measurable effect 
on the geometry of the problem. One might try to interpret a well known model 
for wave propagation in cold plasma, in which the parabolic line of an elliptic- 
hyperbolic system is a parabola with vertex at the origin of [55] , [105] , as a 
fixed gauge of a Hodge system on extended P^. 

Unfortunately, there would be some problems with interpreting the cold plasma 
model in terms of waves on extended P^. First, it does not appear to simplify the 
original problem, although it may illuminate certain similarities in the analysis of 
the two systems; compare [HT], [S2j, and [9S]. Second, the projective group consists 
of non-Euclidean motions, which are noninertial and thus lacking in obvious physical 
meaning, so it is unclear how to interpret the mathematical relation between the 
two systems in terms of a satisfying physical theory. Finally, the gauge groups of 
Sec. 5.1 act "upstairs" on a fiber bundle of states. The symmetry group under which 
the Laplace-Beltrami equations are invariant acts "downstairs" on the underlying 
metric, in the manner of the gauge group of general relativity. Thus it is reasonable 
to look for physical analogies of Sec. 2 among theories in which the relevant bundle 
is soldered to the base space rather than being related to the base space only by a 
projection, as in plasma physics, the Born-Infeld model, and other electromagnetic 
theories. 

Symmetry groups for differential operators on Riemannian-Lorentzian metrics 
other than extended P^ have been computed |69| . In particular, the symmetry 
group for operators of Tricomi type which satisfy a pure power law has been shown 
to correspond to a group of local conformal transformations with respect to the 
underlying metric away from the metric singularity; moreover, the group extends 
across the singular surface on which the metric changes type |102j . 



For the reasons cited in the preceding section, it is reasonable to look to relativity 
for physical analogies of Sec. 2. We therefore include a brief survey of signature 
change in the recent physics literature. 

1. Special relativity: The wave equation on Minkowski space-time, in a reference 
frame rotating with constant angular velocity u with respect to another reference 
frame, is expressible in cylindrical coordinates (p, ip, z) as the elliptic-hyperbolic 
equation [lllj 
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2. Quantum cosmology: These examples arise from the (controversial) Hartle- 
Hawking hypothesis |46| . that the universe might have originated as a manifold 
having Euclidean signature and subsequently undergone a transition to a model 
having Lorentzian signature across a hypersurface which was space-like as seen 
from the Lorentzian side. (Note that certain metrics which are called Euclidean by 
physicists would be called Riemannian by geometers; see footnote 2 of [36].) Some 
2-dimensional variants are [lllj : 

i) continuous change of signature: 

ds"^ -tdt'^ + dz^] 

ii) discontinuous change of signature: 

ds^ = -z^'^dt^ + dz'^\ 
Hi) continuous change of signature with a curvature singularity: 

ds^ = -zdt^ + dz^; 

see also [27]. 

These examples have obvious higher-dimensional analogues. 
Alternatively, the metric might change from Lorentzian to Klcinian signature 
across the line z = : 

ds^ ^ -dt^ + dx^ + dy^ + zd??; 
this 4-dimensional model has been studied in, e.i?., [3j. 

In connection with the distinction between examples ii) and iii)^ we note that 
operators on Riemannian-Lorentzian metrics which degenerate rather than blow up 
at the change of signature have been studied by mathematicians as well as physicists 
[19]; see also [18^ and the discussion in Sec. 5 of [102j . 

It is interesting for mathematicians that although many of the criticisms of Har- 
tle and Hawking's widely discussed proposal concern physical predictions that it 
implies, controversy also arises from mathematical ambiguities in geometric anal- 
ysis on mixed Riemannian-Lorentzian metrics and from the variety of potential 
singularities of such metrics. See [99j and the references therein for a recent dis- 
cussion of the physical predictions; see, e.g., [M], [22], [S], [32], [J^, [50], [51], and 
[55] for discussions of the mathematical ambiguities. It has been observed [3D] that 
singularities similar to those that are associated with the Hartle-Hawking transition 
can also arise in classical relativity, as in Example 1, above; see also [25], [28] . 
[48], [60], and [61]. 

3. Repulsive singularities in 4-dimensional extended supergravity |35j : Mathe- 
matically, this model is essentially a combination of 2ii) and 2iii); 

4. Binary black hole space-times with a helical killing vector 58J: This model is 
a generalization of example 1; see also [52] and [116j . 

5. Brane worlds: A hrane is a submanifold of a hulk, or higher-dimensional space- 
time. Traditionally, branes have been represented as uniformly time-like, but it has 
been observed ^75j that they need not be time-like everywhere and provide a natural 
context for signature change. In fact, mixed Euclidean-Lorentzian branes can be 
constructed in such a way that both the bulk and the brane are regular. However, 
if viewed from within the brane, the change of signature may appear as a curvature 
singularity [36], [76j . For certain choices, this provides an elegant kinematic model 
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for both the apparent big bang singularity and the apparent accelerated expansion 
of the universe; in particular, no hypothesis of dark energy is needed to account for 
accelerated expansion [77j . 

Signature change has also been investigated in the context of spinor cosmology; 
see [120] for a recent example. 

There is a substantial literature on constructing analogies for the dynamical 
equations for light in curved space-time using equations from models of condensed 
matter, including (but by no means limited to) Bose-Einstein condensates with a 
sink or a vortex [7], [66], |121j - [125] . These lead to analogies between acoustic waves 
in matter and wave equations on mixed Riemannian-Lorentzian manifolds which 
are reminiscent of those between gas dynamics and extremal surfaces discussed in 



The idea behind these models is elegantly simple: An acoustic wave has a relation 
to the flow in which it propagates which is analogous in some crucial ways to the 
relation between a light wave and the ambient space-time. This analogy yields a 
kinematic model for certain relativistic effects. For example, if the flow becomes 
supersonic, an acoustic wave emitted downstream from a listener will be trapped 
in an analogous way to the trapping of light inside a black hole with respect to an 
external observer. 

Such analogies can be traced back at least to the flow metrics of [12] and possibly 
to the electrodynamics of [38j ; a review is given in [6J . 

An elliptic-hyperbolic system associated with the Einstein evolution equations 
is studied in 0] . But in that example an elliptic gauge-fixing condition is coupled 
to hyperbolic evolution equations on a Lorentzian metric. Because the signature 
of the metric is fixed, the resulting system is qualitatively different from the cases 
considered here. 

An important consideration for physical applications is whether the elliptic- 
hyperbolic differential operator is of real principal type, in which case the principal 
symbol of the operator is real-valued and no complete null bicharacteristic can be 
trapped over any compact subset of the domain. Because a null bicharacteristic is 
an integral curve of a Hamiltonian system canonically associated to the principal 
symbol, the major analytic properties for operators of real principal type depend 
only on the principal symbol, and not on the form of the lower-order terms; see the 
concluding remarks in Sec. 3 of [lOlj. 




If the operator is of real principal type, ideas from microlocal analysis can be 
applied to construct a natural theory of boundary regularity that is applicable to 
the elliptic-hyperbolic case jlOOj ; see also [39] and [4^. A major difference between 
most of the physical examples listed here and the geometric examples of Sees. 1-5 
is that the latter operators are not of real principal type; as a result, microlocal 
arguments appear to fail and what can or cannot be said about solutions tends to 
depend delicately on the precise form of the lower-order terms. A typical example 
is Theorem 3 of |41j . In particular, this kind of dependence prevents the derivation 
of uniqueness theorems by the expected arguments; see Sec. 5.1 of [98j for an 
illustration. 



Sec. 4.2. 




Acknov^rledgment . I am grateful to Prof. K. R. Payne for acquainting me with 
the accepted terminology for equations of the form (fT4|) . Confusion arises because 
M. V. Keldysh was not the first mathematician to study such equations. 
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